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Summary. We present the results of mesoscopic dissipative particle dynamics 

(DPD) simulations of coupled electrohydrodynamic phenomena on the micro- and 
Q . nanoscale. The effects of electroosmotic flow and slippage combined with polyelec- 

& . trolyte electrophoresis are investigated in detail, taking full account of hydrodynamic 

and electrostatic interactions. Our numerical results are in excellent agreement with 

analytical calculations. 



Ch ' 1 Introduction 

O" 

Microfluidic devices like bio-MEMS (micro-electronical-mechanical-systems) and 
bio-NEMS (nano-electronical-mechanical-systems) have attracted broad interest 
over the last years due to their huge potential in biotechnology. The flow profiles 
^ ' in such micro- or nanosized devices are strongly influenced by the properties of the 

**i , boundaries due to the large surface-to-volume ratio in these systems. Surface char- 

acteristics like the wetting behavior and/or slippage have a dramatic effect on the 
microscopic flow, leading to sometimes unexpected behavior. 
One particularly important mechanism is electroosmotic transport: in contact with 

C* — ' a liquid, many materials commonly used in nanotechnology (e.g., polydimethylsilox- 

^^ \ ane (PDMS)) become charged due to ionizations of surface groups [Isr91]. As a con- 

sequence, surfaces are often covered by a compensating counterion layer [Hun89]. In 
an external electric field, the ions are driven in one direction, dragging the surround- 
ing solvent with them. As a result, a flow is induced in the fluid, the electroosmotic 
flow (EOF). This electrokinetic effect has numerous consequences. For example, it 

|r% ' alters drastically the migration dynamics of mesoscopic objects like polyelectrolytes 

or colloids [VioOO] . In microchannels, the EOF generated at the channel walls induces 
a total net flow, which is technologically attractive because it can be controlled and 
manipulated more easily on the submicrometer scale than pressure- or shear-driven 
flow. 

One important application of microchannels is to separate different fragments of 
biological matter like DNA by their length for sequencing or further manipulation. 
High molecular weight polyelectrolytes cannot be separated by naive electrophoresis 
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in free salt solution due to the fact that the electrophoretic mobility becomes length 
independent for long chains [VioOO]. In many separation methods, the polyelec- 
trolytes are driven through micro- or nanostructured environments to overcome this 
problem, e.g., a disordered gel (in gel electrophoresis), or structured microchannels. 
The presence of boundaries alters the dynamical behavior of the macromolecules 
drastically [MMT07], and their dynamical behavior results from a complex interplay 
of electrostatics, hydrodynamics, and confinement effects on the molecules. 
Our research focuses on the investigation of the explicit coupling of electrohydrody- 
namic effects in the nanometer scale in confined geometries by coarse-grained meso- 
scopic simulations [SSH09, SS10]. We use Dissipative Particle Dynamics (DPD), 
which is a popular mesoscopic method in mesoscopic simulations. The results for 
polyelectrolyte electrophoresis in small microchannels are related to the experiments 
published by Mathe et al. in [MMT07]. In particular, we investigate in detail the 
influence of the electroosmotic flow on the total polyelectrolyte mobility. We find an 
excellent agreement between theory and numerical results. 



2 Dissipative Particle Dynamics 

Dissipative Particle Dynamics (DPD) was originally developed by Hoogerbruggc 
and Koelman [HK92] as a combination of Lattice Gas Automata characteristics and 
Molecular Dynamics methods. Compared to atomistic Molecular Dynamics simula- 
tions, this method gives access to much longer time- and length scales and is there- 
fore suited to study the long-time behaviour of soft matter systems and transport 
phenomena. It is coarse-grained, momentum-conserving, and creates a well-defined 
canonical ensemble. 

The basic DPD equations are given by the forces on one particle, which involve 
two-particle interactions that are given by 



Fr D =]TFg+Fg+F 



i¥=3 



with a conservative force F*j 



») 



Fg = -V l] U i] {r lj ), (2) 



a dissipative force F^ which reads 

Ffj = -fDPDUD(rij)(fij ■ Vij)fij (3) 

with the friction coefficient jdpd and an additional random force F^ yielding 

F^=acj R (r ij )Cijfij. (4) 

where the weighting function is given by 

uo(rij) = [Wfl(rij)] = ujDPD{rij), (5) 

such that a canonical ensemble is sampled at equilibrium [GW97, EW95]. In DPD 
simulations, the conservative forces are often taken to have a certain soft shape 
[GW97]. Here we only use the DPD thermostat as described above. The random 
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number £ has zero mean and unit variance and is symmetric, Qj — Qi, to ensure the 
conservation of momentum and the weighting function is arbitrary and often chosen 
linear [GW97]. It depends on the interparticle distance Tij and the cut-off radius r c . 
The strength of the interaction is steered by the distance of the particles with 

while the amplitude of the Gaussian white noise in Eqn. (4) is given by 

o = 2jDPDkBT (7) 

with the Boltzmann constant fcs and the temperature T. Eqn. (1) can be integrated 
by an ordinary Molecular Dynamics integration scheme like the Velocity- Verlet al- 
gorithm [FS01] which is used in the ESPResSo package. 



3 The software package ESPResSo 

All simulations in this work have been carried out using extensions of the soft- 
ware package ESPResSo (An Extensible Simulation Package for Research on Soft 
matter) [AML05]. ESPResSo was mainly developed for coarse-grained mesoscopic 
simulation approaches. One of the advantages of this program is its high perfor- 
mance MPI-parallelization implemented for simulations on supercomputers. Users 
can furthermore change and extend the program code to adopt it for their own pur- 
poses. ESPResSo incorporates several simulation techniques like Lattice-Boltzmann, 
Dissipative Particle Dynamics, Stochastic Dynamics as well as pure Molecular Dy- 
namics techniques. Another feature are the several implemented electrostatic algo- 
rithms like MMM1D, MMM2D, MMM3D, P3M, ELC and screened Debye-Hiickel 
potentials, which allow to choose between the fastest calculation methods available. 
In summary, ESPResSo provides a well founded basis code for high performance 
computing on parallel clusters. The steering of the simulations is based on a TCL 
(Tool Command Language) script. For starting a simulation, no explicit knowledge 
of implementation details is needed. Even newcomers in the methods of computer 
simulations can succesfully run a simulation after a short time. 
Several tools for analysis are additionally included in the program. ESPResSo is 
under public license and free to download [AML05]. Users that develop new ideas 
are invited to submit their source code, written in the programming language C 
to be incorporated after testing in the newest release version. The development of 
ESPResSo continues and a number of extensions is already planned. 
We ran our simulations on the NEC SX-8 Cluster at the High Performance Com- 
puting Center Stuttgart. The number of computing nodes has normally been chosen 
to 4 which corresponds to 32 CPUs. Each job has got a typical runtime of 4 hours 
by using the concept of setting checkpoints. 



4 Polyelectrolyte electrophoresis in microchannels 

4.1 General theory 

We consider for simplicity a planar slit channel with identical walls at z = doL/2, 
exposed to an external electric field E x in the x direction. The electrostatic potential 
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<t>{x,y,z) then takes the general form <P(x,y,z) = ip(z) + E x x + const, where we 
can set tp(Q) = for simplicity. The electrolyte in the channel is taken to contain n 
different ion species i with local number density pt(z) and valency Zi, which results 
in a net charge density p(z) = ^"_ (Zie)pi(z). The electric field then generates a 
force density f x (z) — p(z)E x in the fluid. Comparing the Poisson equation for the 
electrostatic potential rp, 

d\{z) _ P {z) 

dz 2 e r 

(where e r is the dielectric constant), with the Stokes equation 



(«) 



Vs^^- = -Mz) = -p(z)E a (9) 

(with the shear viscosity rj 3 ), one finds immediately d zz v x (z) = d zz ip(z) (e r E x /r/ s )- 
For symmetry reasons, the profiles v x and if) must satisfy the boundary condition 
d z v x \ z= o = d z tp\ z= o — at the center of the channel. This gives the relation 

e r E 
v x (z) = ^— —ip(z)+v EOF , (10) 

where we have used ip{jS) =0 and identified the fluid velocity at the center of the 

channel with the EOF velocity, v x (0) = f EOF . We further define ipB ■— V>(±zb) (for 

no-slip boundaries, ipB is the so-called Zeta-Potential [Hun89, SS10]). 

The roughness of the channel boundaries is included in the partial slip boundary 

condition. 

5b d z v(x)\ ZB — v x (z)\ ZB , (11) 

where v x (z) denotes the component of the velocity in x-direction evaluated at the 
position zb of the so-called "hydrodynamic boundary" . This boundary condition is 
characterized by two effective parameters, namely (i) the slip length 8b and (ii) the 
hydrodynamic boundary zb- We note that the latter does not necessarily coincide 
with the physical boundary. Inserting the partial-slip boundary condition for the 
flow, Eqn. (11), we finally obtain the following simple expression for the electroos- 
motic mobility, 

Meof = v EOF /E x = p° EOF (1 + k 5b), (12) 

where we have defined the inverse 'surface screening length' 

d z ip 
k--=T-—\z=±z b , (13) 

and A«eof la ^ ne well-known Smoluchowski result [Hun89] for the electroosmotic 
mobility at sticky walls, 

Cf = -£rip B /Vs- (14) 

The remaining task is to determine the screening parameter k. If the surface 
charges are very small and the ions in the liquid are uncorrelated, it can be calculated 
analytically within the linearized Debye-Hiickel theory [Isr91]. The Debye-Hiickel 
equation for the evolution of the potential ijj in an electrolyte solution reads d zz ip — 
K 2 D ip with the inverse Debye-Hiickel screening length 



KD = \ e r k«T ■ (15) 



Mesoscopic simulations of polyelectrolyte electrophoresis in nanochannels 5 

where p^o is the density of ions i far from the surface. It is solved by an exponentially 
decaying function, 

i>(z)<x(e KDZ + e~ KDZ -2). (16) 

Inserting that in Eq. (13), one finds k = Kd, i.e., the surface screening length is 
identical with the Debye screening length. 

Unfortunately, the range of validity of the Debye-Hiickel theory is limited, it 
breaks down already for moderate surface potentials i\>b and/or for highly concen- 
trated ion solutions. Nevertheless, the exponential behavior often persists even in 
systems where the Debye-Hiickel approximation is not valid. For high ion concen- 
trations a Debye-Hiickel-type approximation can still be used in a wide parame- 
ter range, if kd is replaced by a modified effective screening length [BKN05]. For 
high surface charges, analytical solutions are again available in the so-called 'strong 
coupling limit', where the profiles are predicted to decay exponentially with the 
Gouy-Chapman length [MN01]. This limit is very special and rarely encountered. 
At intermediate coupling regimes, the decay length must be obtained empirically, 
e.g., by fitting the charge distribution p(z) to an exponential behavior, which is 
characterized by the same exponential behavior than tp(z) by virtue of the Poisson 
equation, 

^(Z ie )^)oc^^oc(e"'+e-"'). (17) 

i = q 

Assuming that the electrophoretic velocity of the polyelectrolyte 

v p (x) = fi e E x (18) 

is influenced by the electroosmotic mobility of the counterions /j,eof, a total mo- 
bility has to be defined /i t = /i EOF + fi e which describes the overall mobility of the 
electrophoretic object. Putting everything together, the total net electrophoretic 
mobility fM of a polyelectrolyte in the channel can be expressed in terms of the 
electroosmotic mobility /i EOF as 

Meof MeofU + ^b) 

where the ratio /n e /At EOF depends only weakly on the ionic strength of the electrolyte 
and the slip length of the surface. The main effect of slippage is incorporated in the 
factor (I + kSb)' 1 ! [SS10]. 

4.2 Simulation details 

We have studied the electrophoresis of a single charged polymer of length N = 20 
in electrolyte solutions, confined by a planar slit channel with charged walls. All 
particles, polymer, solvent and ions, are modeled explicitly. We use a simulation 
box of size (12a x 12a x 10a) which is periodic in x- and {/-direction and confined by 
impermeable walls in the 2-direction. The walls repel the particles via a soft repulsive 
WCA potential of range a and amplitude e. (Hence the accessible channel width for 
the particles is actually L z — 8a). Ions and monomers repel each other with the same 
WCA potential. In addition, chain monomers are connected by harmonic springs 

1 2 

U harmonic — TT v *J '"O/ \ "/ 
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with the spring constant k — 25e/a 2 and ro = 1.0a. Neutral solvent particles have 
no conservative interactions except with the walls. The wall contains immobilized, 
negatively charged particles at random positions. Every second monomer on the 
polyelectrolyte with 20 beads carries a negative charge resulting in 10 charged 
monomers. The solvent contains the positive counterions for the walls and the 
polyelectrolyte, and additional (positive and negative) salt ions. All charges are 
monovalent, and the system as a whole is electroneutral. In addition to their other 
interactions, charged particles interact via a Coulomb potential with the Bjer- 
rum length Xb = e 2 /47re r fcsT = l.Ocr, and they are exposed to an external field 
E x = — l.Oe/ea. Specifically, we have studied systems with a surface charge den- 
sity of o a — — 0.208ea~ 2 which corresponds to 30 charged particles per wall. In 
a recent publication [SSH09], we have shown that this corresponds to the 'weak- 
coupling regime', i.e., the regime where the Poisson-Boltzmann theory is valid. The 
electrostatic coupling constant [MN01] S — 2ir Z 3 ; A s a a (Z — 1 is the valency of the 
cations), which gives the strength of electrostatic interactions between the surface 
and the ions compared to thermal energy, is close to unity, S ~ 1.3. In this earlier 
work, we have also studied the effect of using homogeneously charged wall instead 
of discrete embedded charge, and found the differences to be negligible [SSH09]. 
The use of discrete embedded charges has practical advantages in our simulation 
code, which is why we use them here. The total counterion density was p = 0.06a" -3 
and the salt density varied between p s =0.05625, 0.0375, 0.03, 0.025, and 0.015a -3 . 
In molar units, this corresponds to 0.272, 0.181, 0.145, 0.121 and 0.072 mol/1, if 
we identify As ~ 0.7 nm, i.e., the Bjerrum length in water at room temperature 
[VioOO]. 

We use DPD simulations with a friction coefficient 7dpd = 5.0a _1 (me) 1 ' 2 . (More 
precisely, we only use the DPD thermostat, not the soft conservative DPD forces. All 
conservative forces in our system are either WCA forces, spring forces, or Coulomb 
forces as explained above). The density of the solvent particles was p = 3.75a -3 , 
and the temperature of the system was T — l.Oe/fcs. For these parameters, the 
shear viscosity of the DPD fluid - as determined by fitting the amplitude of Plane 
Poiseuille flows [SAF08] - is given by rj s = (1.334 ± 0.003)a" 2 (me) 1/2 . The DPD 
timestep was St = 0.01a(m/e) 1 ^ 2 . 

Tunable-slip boundary conditions were used with friction coefficients 71, =0.1, 
0.25, 0.5, 0.75, 1.0, and 6.1a~ 1 (me) 1 ' 2 . The range of the viscous layer was z c = 2.0a. 
Only the solvent particles interact with the tunable-slip boundaries. By performing 
Plane Poiseuille and Plane Couette flow simulations with the above given param- 
eters, the slip length 8b and the hydrodynamic boundary positions zb can be de- 
termined independently [SAF08]. The hydrodynamic boundary position is found at 
\zb\ — (3.866 ± 0.266)a in all simulations. 

4.3 Results 

Polyelectrolytes in absence of external fields 

To understand the behaviour of polyelectrolytes in salty solution we first studied 
the dynamics in absence of external fields. The general theory [VioOO] implies that 
hydrodynamic interactions between the monomers are important to describe the 
dynamic behaviour adequately. This results in the so called Zimm dynamics [VioOO] . 
If external electric fields are present, hydrodynamic interactions between the charged 
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monomers are screened due to the presence of salt. The two rivaling effects which 
are the electrophoresis of the polyelectrolyte as well as the electroosmosis of the 
mobile counterions in opposite direction lead to a crucial screening of hydrodynamic 
interactions which is described by Rouse dynamics [VioOO] . 

A powerful tool to investigate the underlying dynamics in experiments as well as in 
computer simulations is given by the dynamic version of the structure factor which 
is defined by 

S ( k '^ = Jf S < e lk(R > (t) - R ' (to)) > (21) 

i,i 
with the actual monomer position R; respectively Rj. For the inverse length scale 
1/Rg <C k <C 1/ao with the gyration radius R g [DE86] and the smallest microscopic 
length ao and the finite time interval if, <C t <C r after the ballistic time tb and the 
longest relaxation time r, the dynamic structure factor obeys the following scaling 
relation [DK93] 

S(k,t) = S{k,0)f(k z t) (22) 

which depends on the the parameter z. This parameter differentiates between the 
different regimes. For the Zimm dynamics it is given by 



and for Rouse dynamics 



z = 3 (Zimm-Regime) (23) 



2 = 2 + l/i/ (Rouse-Regime) (24) 



with the Flory paramter v which describes the inverse fractal dimension of the 
polyelectrolyte [DE86] . 

We have compared an uncharged polymer and a half charged polyelectrolyte in a 
salt solution of concentration p s = 0.05a -3 and a solvent density of p — 3. Oct -3 with 
free periodic boundary conditions in all dimensions. All the other parameters are in 
agreement to the above described simulation details. Fig. 1 presents our simulation 
results for an uncharged chain and the polyelectrolyte. Our results indicate that 
although external fields are absent, polyectrolytes can be described best by Rouse 
dynamics whereas uncharged polymers can be described by Zimm dynamics. More 
results on that topic can be found in [Smi09, SS08]. 

Polyelectrolyte electrophoresis in microchannels 

In this section our numerical results are presented which indicate the drastic influ- 
ence of the electroosmotic flow on the total mobility of the polyelectrolyte. Larger 
slip lengths even enhance this effect. 

Fig. 2 presents the results of the ionic distribution for a salt concentration of 
p s = 0.05625a -3 and a counterion density of p — 0.0525a -3 with a surface ion den- 
sity of a s = 0.208a -2 in presence of the polyelectrolyte. Due to the large number of 
cations in comparison to the anions, the cation density in the middle of the channel 
is increased (inset of Fig. (2)). Calculating the ionic difference Ap — p c — p a as the 
difference of the cationic and anionic density yields 

2 N 

£ si) ( z) = -J2^Pi(z) = -^A P (25) 
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Fig. 1. Dynamic structure factor S(k,t) for an uncharged chain (red) in correspon- 
dence to a half charged chain (blue) for 0.2 < ka < 0.3 and salt concentration 
p s = 0.05er -3 . The timescale for the uncharged chain is < t/r < 180 in con- 
trast to < t/r < 300 for the half charged polyelectrolyte. Both chains consist of 
N — 50 monomers. Left side: Zimm-scaling with z = 3. Right side: Rouse-scaling 
with z — 3.5. 



Q. 

< 




Fig. 2. Distribution of the ionic difference Ap = p c — p a for an exemplary salt 
concentration of p B — 0.05625a -3 with a surface ion density of a s = 0.208cr -2 . 
The red line corresponds to the exponential fit of Eqn. (26) with an effective inverse 
screening length of K — 2.305±0.025ct _1 . Inset: Distribution of cations (salt cations 
and counterions) with green symbols and anions (salt anions) with black symbols 
for an the exemplary salt concentration for the above given parameters. 
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Table 1. Fitted inverse screening lengths k for the different salt concentrations p s 
and a fixed counterion density of p = 0. 0525a _ . 



Ps\v- 3 ] 


nia- 1 ] 


0.015 


1.996 ± 0.041 


0.0225 


2.011 ± 0.049 


0.03 


1.983 ±0.041 


0.0375 


2.182 ± 0.047 


0.05625 2.305 ± 0.025 



with the unique valency \Z\ [SS10]. As we have mentioned in section 4, the applica- 
tion of the Debye-Hiickel theory for non electroneutral bulk systems is not valid and 
has to be replaced by an effective inverse screening length n. Thus we have used the 
following fit function 

Ap = Ap () (e- KZ + e KZ ) + c (26) 

to determine the effective inverse screening length k, with the fit parameters Apo, c 
and kd- The red line in Fig. 2 shows that the fit describes the ionic difference ad- 
equately. The fitted parameter values of k for the different salt concentrations are 
presented in Tab.l. Significant differences between the inverse screening lengths are 
only observed for salt concentrations p s > 0.0375<j~ 3 . Thus, for lower salt concen- 
trations, the screening of electrostatic interactions is mainly effectuated by the much 
larger number of counterions instead of the salt ions [SS10]. 

The corresponding electroosmotic flow profiles for an exemplary salt concentra- 
tion of p s — 0.05625ct~ 3 are shown in Fig. 3. The different magnitudes correspond 
to varying slip lengths. Larger slip lengths enhance the flow profile drastically. All 
points are in good agreement to the integrated analytical expression of Eqn. (26) 
in terms of the Stokes equation (9) with partial-slip boundary conditions. Thus the 
description of the electroosmotic flow in terms of the Stokes theory is valid. 

Fig. 4 compiles our numerical results for the electroosmotic mobility for all salt 
concentrations and slip lengths. They are in very good agreement to the theoreti- 
cal prediction of Eqn. (12), where p g EOF has been determined independently by a 
linear regression for each salt concentration. It is worth noting that the presence of 
the polyelectrolyte does not perturb the amplitude of the electroosmotic flow. 
After investigating the electroosmotic flow of the solvent, we focus on the dynam- 
ics of the polyelectrolyte. The monomer distribution for a salt concentration of 
p s — 0.05625<7~ 3 is presented in Fig. 5. The fit shows that the distribution is dom- 
inated by a peak in the middle of the channel with a variance of Var ~ 2.28a. 
Thus the main drag of the electroosmotic flow on the polyelectrolyte is exerted in 
the middle of the channel, and the rapid variations of the electroosmotic flow in 
close vicinity to the boundaries (Fig. 3) have little influence on the polyelectrolyte 
mobility. Estimating the total mobility by assuming a constant plug-like flow profile 
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Fig. 3. Exemplary flow profiles for a salt concentration p s = 0.05625a -3 for varying 
slip lengths (from bottom to top: S B = (14.98, 5.66, 2.63, 1.77, 1.29, 0.00)a.) The 
black lines are the integrated theoretical predictions in terms of the Stokes equation 
(9) with a fitted inverse screening length of k = 2.305a -1 . 



for the fluid in the middle of the channel from z — —1.52a to Z — 1.52a due to the 
variance of the monomer distribution therefore seems reasonable. 
The influence of the electroosmotic flow on the polyelectrolyte can be investigated 
by regarding static properties like the radius of gyration R^ = (l/2iV 2 ) ^\ ._ < 

(Ri - Rj) 2 > and the end-to-end radius with R\ =< (R N - Ri) 2 > [DE86]. The 



Table 2. Radius of gyration R g and end to end radius R e for a polyelectrolyte with 
iV = 20 monomers for different salt concentrations p a . 



Ps[o~ 3 ] 


R g [a] RM 


0.015 


3.2218 ± 0.047 10.6480 ± 0.0314 


0.0225 


3.1661 ± 0.0041 10.2736 ± 0.0266 


0.03 


3.1486 ± 0.0451 10.1777 ± 0.0292 


0.0375 


3.1279 ± 0.0045 10.0819 ± 0.0287 


0.05625 3.0825 ± 0.0045 9.8331 ± 0.0280 
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Fig. 4. Ratio /U EOF /jU EOF plotted against Sb K for the different salt concentrations 
and screening lenghts given in Tab. 1 with counterions of density p = 0.0525cr . 
The blue line is the theoretical prediction of Eqn. (12) with slope 1 + SbK. 



values of these parameters are shown in Tab. 2. Both properties decrease with larger 
salt concentration due to a more pronounced screening of electrostatic interactions 
which is in accordance to standard theories [VioOO]. Although the values for the 
end to end radius are quite large, in all cases the box dimensions in the unconfined 
directions with b x ,y — 12cr are larger than the maximal average extension of the 
chain. Possible explanations for these large ratios between the end-to-end and gy- 
ration radii could be shear-induced elongation [BKN05] or direct squeezing due to 
the presence of the channel walls. 

The total mobility of the polyelectrolyte for varying boundary conditions is finally 
presented in Fig. 6. The theoretical prediction of Eqn. (19) is in good agreement 
with the numerical results with the fitted ratio fj ie /fjP EOF = — 3.778±0.128. With no- 
slip boundaries (Sb ~ 0) one obtains ordinary behaviour where the polyelectrolyte 
migrates in opposite direction to the electroosmotic flow. In the presence of slip, 
the absolute mobility may become negative if the amplitude of the electroosmotic 
flow exceeds a critical value given by a combination of the inverse screening and slip 
lengths, and if the sign of the wall charges is identical to the net charge of the poly- 
electrolyte. If the wall is oppositely charged to the polyelectrolyte, slippage effects 
should even enhance the total velocity of the polyelectrolyte. 

The right inset of Fig. 6 shows the total displacement of the chains center of mass 
for various slip lengths and a salt concentration of p s = 0.05625<r -3 . In nearly all 
cases except for 8b ~ 0, the total mobility of the polyelectrolyte is negative with 
\fi p \ -C IMbof I' indicating negative values of fit/ 'n EOF . This can only be explained 
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Fig. 5. Monomer appearance probability inside the channel for a salt concentration 
of p s = 0.05625a- 3 . 



in terms of the drag force of the electroosmotic flow and the surplus of cations in 
close vicinity to the boundaries (Fig. 2) and is in agreement to the results derived 
above. 

To summarize, the assumptions and calculations of section 4 are in good agree- 
ment to the presented numerical results in this section. The total mobility of the 
polyelcctrolytc can therefore be adequately described by Eqn. (19). 



5 Summary 

We have presented mesoscopic DPD simulations of polyelcctrolyte electrophoresis 
in narrow microchannels, taking full account of hydrodynamic and electrostatic in- 
teractions. We have shown that the product of the inverse screening length k and 
the slip length 8b massively influences the electroosmotic flow and therefore the 
total mobility of the polyelectrolyte. Thus the characteristics of the boundaries have 
to be taken into account for a proper description of the polyelectrolyte migration 
dynamics. For certain parameter sets, even a negative mobility can be achieved. All 
our numerical results are in good agreement to the analytical derived results. 
In summary, only a combination of electroosmotic, electrophoretic, electrostatic and 
slippage effects does describe the total mobility of polyelectrolytes in microchannels 
adequately. Our simulations indicate and explain total negative mobilities due to 
boundary effects which also have been observed in recent experiments [MMT07]. 
The characteristics of the channel walls could be used to significantly enhance flow 
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Fig. 6. Ratio im/\h e 



plotted against 8b k for all salt concentrations. The black 



line is the theoretical prediction of Eqn. (19) with absolut values of |/Lt BOF |. In the 
limit 5bk — >• oo, the total mobility of the polyelectrolyte is equal to the electroos- 
motic flow. The ratio fj, p /fj? EOF has been fitted to —3.778 ± 0.128. Negative values 
of /xt/|/i EOF | indicate absolute negative total mobilities of the polyelectrolyte. In- 
set: Total displacement of the polyelectrolytes center of mass for different boundary 
conditions and a salt concentration of p s — 0.05625ct~ 3 . The total mobility becomes 
negative if the relation |/i e | ^ IMeof I * s fulfilled. The lines correspond from top to 
bottom to the slip lengths S B » (0.00, 1.292, 1.765, 2.626, 5.664, 14.98)cr. Thus larger 
slip lengths indirectly enhance the total mobilty of the polyelectrolyte. 



profiles which offers the possibility to reduce the time needed for polymer migration 
or separation techniques. This could be an important aspect for future applications 
in microchannels or micropumps to accelarate the measuring time in experiments. 
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